Model 4

Description of the model


This fourth model is a new variation. It still compute the mass-action on the numerator, but attempts to make the flux saturate with the part on the denominator with the part of \(1 + {...}\), where we use handling time and the whole biomass consumption of the predator. \[\begin{align} F_{ij}^{real} = \frac{F_{ij}^{theoretical}}{1+h_j * \sum_{i=1} ^{i}{F_{j}^{theoretical}}} && F_{ij}^{theoretical} &= \alpha_{j} * B_i * \frac{B_j}{M_j} \\ \end{align}\]

où
- \(\alpha\) est un paramètre spécifique à un prédateur (mais le prédateur peut se retrouver dans plusieurs réseaux différents)
- \(B_i\) est la biomasse de la proie dans un réseau spéficique
- \(B_j\) et \(M_j\) sont la biomasse et la bodymass d’un prédateur spéficique. La biomasse du prédateur est la même dans un réseau, mais peut varier d’un réseau à l,autre, alors que son bodymass sera le même pour chacun des réseaux.
- \(h_j\) est une variable déjà calculée.

Summary table

Summary table model 4
mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
a_pop -6.971942 0.0060072 0.3347252 -7.610802 -7.197355 -6.979476 -6.752147 -6.294061 3104.788 1.0000643
a_sd 2.714704 0.0049257 0.2284576 2.296760 2.556516 2.702836 2.857735 3.204793 2151.157 1.0016759
sigma 1.660158 0.0003105 0.0322949 1.598286 1.637839 1.659591 1.681981 1.723494 10815.021 0.9997148

Exploring parameters posterior distribution

Respective predator alphas

Respective predator handling times

One-one plot of simulation against data